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Abstract 



The Cryogenic Dark Matter Search (CDMS) employs Ge and Si detectors to 
search for WIMPs via their elastic-scattering interactions with nuclei while 
discriminating against interactions of background particles. CDMS data give 
limits on the spin-independent WIMP-nucleon elastic-scattering cross-section 
that exclude unexplored parameter space above 10 GeV c~ 2 WIMP mass and, 
at > 84% CL, the entire 3<r allowed region for the WIMP signal reported by 



1 



the DAMA experiment. 

PACS numbers: 95.35.+d, 14.80.-j, 14.80.Ly 
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Extensive evidence indicates that a large fraction of the matter in the universe is non- 
luminous, nonbaryonic and "cold" - nonrelativistic at the time matter began to dominate 
the energy density of the universe Weakly Interacting Massive Particles (WIMPs) 

are an excellent candidate for nonbaryonic, cold dark matter ||,|2|]. Minimal supersymmetry 
provides a natural WIMP candidate in the form of the lightest superpartner, with a typi- 
cal mass M ~ 100 GeV c~ 2 |5]-§|. WIMPs are expected to have collapsed into a roughly 
isothermal, spherical halo within which the visible portion of our galaxy resides. WIMPs 
scatter off nuclei via the weak interaction, potentially allowing their direct detection P,|10 



The expected spectrum of recoil energies (energy given to the recoiling nucleus during the 



interaction) is exponential with a characteristic energy of a few to tens of keV |jTT|. The 
expected event rate is mo del- dependent but is generically 1 kg -1 d _1 or lower |L0 . 



This Letter reports new exclusion limits on the spin-independent WIMP-nucleon elastic- 
scattering cross-section by the Cryogenic Dark Matter Search (CDMS). The rate of rare 
WIMP-nucleon interactions is constrained by extended exposure of detectors that discrim- 
inate WIMP-induced nuclear recoils from electron recoils caused by interactions of back- 
ground particles [ |T2| , |i"3"| . 

The ionization yield Y (the ratio of ionization production to recoil energy in a semi- 
conductor) of a particle interaction differs greatly for nuclear and electron recoils. CDMS 
detectors measure phonon and electron-hole pair production to determine recoil energy and 
ionization yield for each event. The data discussed here were obtained with two types of de- 
tectors, Berkeley Large Ionization- and Phonon-mediated (BLIP) and Z-sensitive Ionization- 
and Phonon-mediated (ZIP) detectors fl2|- |l7f| . For both types, the drift field for the ion- 
ization measurement is supplied by radially-segmented electrodes on the faces of the disk- 



shaped crystals |T8| . In BLIP detectors, phonon production is determined from the detector's 
calorimetric temperature change. In ZIP detectors, athermal phonons are collected to de- 
termine phonon production and xy-position. Detector performance is discussed in detail 
elsewhere (TJ,[I|Jr^. 

Photons cause most bulk electron recoils, while low-energy electrons incident on the 
detector surfaces cause low-F electron recoils in a thin surface layer ("surface events"). 
Neutron, photon, and electron sources are used to determine efficiencies for discrimination 
between nuclear recoils and bulk or surface electron recoils. Above 10 keV, CDMS detectors 
reject bulk electron recoils with > 99% efficiency and surface events with > 95% efficiency. 
CDMS detectors that sense athermal phonons provide further surface-event rejection based 
on the differing phonon pulse shapes of bulk and surface events [I7|. This phonon-based 
surface-event rejection alone is > 99.7% efficient above 20 keV. 

The 1-cm-thick, 7-cm-diameter detectors are stacked 3 mm apart with no intervening 
material. This close packing enables the annular outer ionization electrodes to shield the 
disk-shaped inner electrodes from low-energy electron sources on surrounding surfaces. The 
probability that a surface event will multiply scatter is also increased. 

The low rate of WIMP interactions necessitates operation at a site with low background- 
particle flux. CDMS detectors are operated beneath 16 meters-water-equivalent overburden, 
which stops the hadronic component of cosmic-ray air showers and reduces the muonic 
component by a factor of 5. A custom, radiopure extension to a modified Oxford S-400 



dilution refrigerator provides a low-background 20 mK volume [19 



Several layers of shielding surround the cryostat. Outermost is a > 99.9% efficient 
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plastic-scintillator veto to detect muons and thus allow rejection of muon-coincident parti- 
cles. Inside the veto, a 15-cm-thick lead shield reduces the background photon flux by a 
factor of 1000. A 1-cm-thick shield made of ancient lead provides additional photon shield- 
ing inside the cryostat p0| . Samples of all construction materials were screened to ensure 



low radioactive contamination. The measured event rate below 100 keV due to photons is 
roughly 60 keV -1 kg -1 d _1 overall and 2 keV -1 kg -1 d _1 anticoincident with veto. 

Neutrons with energies capable of producing keV nuclear recoils are produced by muons 
interacting inside and outside the veto ("internal" and "external" neutrons, respectively). 
The dominant, low-energy (< 50 MeV) component of these neutrons is moderated by a 



25-cm thickness of polyethylene between the outer lead shield and cryostat [21]. However, 



high-energy external neutrons may punch through the moderator. A simulation of these 
neutrons assumes the production spectrum given in 53] and propagates them through the 
shield to the detectors. The accuracy of the simulation's propagation of neutrons is con- 
firmed by the excellent agreement of the simulated and observed recoil-energy spectra due 
to veto-coincident and calibration-source neutrons. A large fraction of the external neutrons 
are vetoed: ~40% due to neutron-scintillator interactions and an unknown fraction due to 
associated hadronic showers. This unknown fraction, combined with a factor of ~4 uncer- 
tainty in the production rate, makes it difficult to accurately predict the absolute flux of 
unvetoed external neutrons. However, normalization-independent predictions of the simu- 
lation, such as relative rates of single scatters and multiple scatters, relative rates in Si and 
Ge detectors, and the shapes of nuclear-recoil spectra, are insensitive to reasonable changes 
in the neutron spectrum. 

Two data sets are used in this analysis: one consisting of 33 live days taken with a 100 
g Si ZIP detector between April and July, 1998, and another taken later with Ge BLIP 
detectors. The Si run yields a 1.6 kg d exposure after cuts. The total low-energy electron 
surface-event rate is 60 kg -1 d _1 between 20 and 100 keV. Four nuclear recoils are observed 
in the Si data set. Based on a separate electron calibration, the upper limit on the expected 
number of unrejected surface events is 0.26 events (90% CL). These nuclear recoils also 
cannot be due to WIMPs. Whether their interactions with target nuclei are dominated 
by spin-independent or spin-dependent couplings, WIMPs yielding the observed Si nuclear- 
recoil rate would cause an unacceptably high number of nuclear recoils in the Ge data set 



discussed below. Therefore, the Si data set, whose analysis is described elsewhere p3 
measures the unvetoed neutron background. 

Between November, 1998, and September, 1999, 96 live days of data were obtained 
using 3 of 4 165 g Ge BLIP detectors. One detector is discarded because it displays a high 
rate of veto-anticoincident low-energy electron surface events, 230 kg -1 d _1 as compared 
to 50 kg -1 d _1 for the other detectors (10 to 100 keV). This detector suffered additional 
processing steps that may have contaminated its surface and damaged its electrodes. Data- 
quality, nuclear-recoil acceptance, and veto-anticoincidence cuts reduce the exposure (mass 
x time) by 45%. To take advantage of close packing, analysis is restricted to events fully 
contained in the inner electrodes, reducing the exposure further by a factor of 2.5 to yield a 



final Ge exposure of 10.6 kg d [24 



Figure |I| shows a plot of ionization yield vs. recoil energy for the Ge data set. Bulk 
electron recoils lie at ionization yield Y ~ 1. Low-energy electron events form a distinct 
band at Y ~ 0.75, leaking into the nuclear-recoil acceptance region below 10 keV. 
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Figure || displays the recoil-energy spectrum of unvetoed nuclear recoils for the Ge data 
set. Only single scatters (events triggering a single detector) are shown; the WIMP multiple- 
scatter rate is negligible. An analysis threshold of 10 keV, well above trigger thresholds, is 
imposed. This choice reduces the data set's sensitivity but simplifies analysis by rendering 
low-energy electron misidentification negligible. The nuclear-recoil efficiency is determined 
using calibration-source neutrons and its stability is monitored using veto-coincident neu- 
trons. Thirteen unvetoed nuclear recoils are observed in the 10.6 kg d exposure between 10 
and 100 keV; this rate is similar to that expected for the WIMP signal claimed by the DAMA 
experiment [p5,2"E|. However, much evidence indicates that the CDMS nuclear recoils are 
caused by neutrons rather than WIMPs. 

Figure [3] displays a scatter plot of ionization yields for multiple scatters. The observation 
of 4 Ge multiple-scatter nuclear recoils is the primary evidence for the neutron interpretation. 
It is highly unlikely that these events are misidentified low-energy electron events. Figures 
|T] and |3] demonstrate excellent separation of low-energy electron events from nuclear recoils. 
Analysis using events due to electrons emitted by the contaminated detector yields an upper 
limit of 0.05 misidentified multiple-scatter low-energy electron events (90% CL). 

All other pieces of evidence are also consistent with the neutron interpretation. First, 
the 4 nuclear recoils observed in the Si data set cannot be interpreted as WIMPs or surface 
events. Second, there is reasonable agreement between predictions from the Monte Carlo 
simulation and the relative observed numbers of 4 Ge multiple scatters, 4 Si single scatters, 
and 13 Ge single scatters. Normalizing the simulation by the 17 total Ge nuclear-recoil events 
yields 2.7 expected Si single scatters and 1.3 expected Ge multiple scatters. Statistically, the 
expected neutron background should result in a less likely combination of Ge single scatters, 
Ge multiple scatters, and Si single scatters 23% of the time. Finally, a Kolmogorov-Smirnov 
test indicates that the deviation between the observed and simulated nuclear-recoil spectral 
shapes would be larger in 28% of experiments. 

The 90% CL excluded region for the WIMP mass M and the spin-independent WIMP- 
nucleon elastic-scattering cross-section a is derived using an extension of the approach of 
Feldman and Cousins ||29|| . The above arguments require accounting for the component of 
the observed Ge single scatters (with energies E[) that is due to the unvetoed neutron flux n. 
This flux is constrained by the number iV m of multiple scatters in Ge and the number of 
nuclear recoils in Si. To determine the 90% CL excluded region in the plane of M and a alone, 
the parameter n is projected out. For a grid of physically allowed values of M, a, and n, the 
expected distribution of the likelihood ratio R = C(Ei, N m , N Si \ a, M,n)/C is calculated by 
Monte Carlo simulation in order to determine the critical parameter R 9Q such that 90% of the 
simulated experiments have R > R$o. Here n is the value of n that maximizes the likelihood 
C for the given parameters M and o and the observations. £ is the maximum of the 
likelihood for any physically allowed set of parameters. The WIMP-nucleon cross-section a 
is converted to a WIMP-nucleus cross-section assuming A 2 scaling with target nuclear mass. 
This scaling is valid for models of supersymmetric WIMPs currently favored The 
90% CL region excluded by the observed data set, with ratio -Rdata, consists of all parameter 
space for which -Rdata < -^90- Figure |] displays the lower envelope of points excluded for all 
values of n. Because all the nuclear recoils may be neutron scatters, a = is not excluded. 

This limit excludes new parameter space for WIMPs with M > 10 GeV c -2 , some of 
which is allowed by supersymmetry [|7J. The data are compatible with the DAMA/NaI-0 ex- 
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elusion limit based on pulse-shape analysis [p7fl. However, these data exclude, at > 84% CL, 



the entire region allowed at 3cr by the DAMA/NaI-1 to 4 annual modulation signal p5 



This region, given by the Vq = 220 km s _1 curve in Figure 4a of Ref. P5| |, is used because 
it is determined solely from the annual modulation signal. The data presented here also 
exclude the analogous 2a allowed region for DAMA/NaI-1 to 2 at > 96% CL |26|. Although 
without theoretical support, non-A 2 scaling may allow the two results to be compatible. 
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FIGURES 




Recoil Energy [keV] 

FIG. 1. Ionization yield (Y) vs. recoil energy for veto-anticoincident single scatters contained 
in the inner electrodes of the 3 uncontaminated Ge detectors. Solid curve: expected position 
of nuclear recoils. Dashed curves: nominal 90% nuclear-recoil acceptance region. Dashed line: 
10 keV analysis threshold. Dashed-dotted curve: threshold for separation of ionization signal from 
amplifier noise. Circled points: nuclear recoils. The presence of 3 events just above the acceptance 
region is compatible with 90% acceptance. 
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FIG. 2. Solid: histogram of nuclear recoils observed in the inner electrodes of the 3 uncontami- 
nated Ge detectors (left-hand scale). Shaded: 10 keV analysis threshold. Dashed: peak-normalized 
nuclear-recoil efficiency (right-hand scale). 
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FIG. 3. Scatter plot of ionization yields for multiple scatters in the 3 uncontaminated Ge 
detectors with at least 1 inner-electrode scatter and both scatters between 10 and 100 keV. Circled 
events are tagged as nuclear recoils in both detectors. Bulk recoils and surface events lie at Y ~ 1 
and Y ~ 0.75, respectively. 
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FIG. 4. Spin-independent a vs. M. The regions above the curves are excluded at 90% CL. 
Solid dark curve: limit from this analysis. Dotted curve: CDMS expected sensitivity (median 
simulated limit) given the observed neutron background. Because the number of multiple scatters 
observed is larger than expected, the limit from this analysis is lower than the median simulated 
limit. Solid light curve: DAMA limit using pulse-shape analysis [27]. Shaded region: DAMA 3<r 
allowed region as described in text pB| . Dashed-dotted curve: Ge diode limit, dominated by [28]. 
All curves are normalized following [o], using the Helm spin-independent form-factor, A 2 scaling, 
WIMP characteristic velocity vq = 220 km s , mean Earth velocity ve = 232 km s , and 
p = 0.3 GeV c~ 2 cm" 3 . 
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